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Abstract 

We present a detailed analysis together with exact numerical calculations 
on one-loop contributions to the neutron electric dipole moment from super- 
symmetry without R parity, focusing on the gluino, chargino, and neutralino 
contributions. Apart from the neglected family mixing among quarks, com- 
plete formulae are given for the various contributions through the quark dipole 
operators, to which the present study is restricted. We discuss the structure 
and main features of the R-parity violating contributions and the interplay 
between the R-parity conserving and violating parameters. In particular, the 
parameter combination fj,*X in , under the optimal parametrization adopted, is 
shown to be solely responsible for the R-parity violating contributions in the 
supersymmetric loop diagrams. While ^*X iu could bear a complex phase, the 
latter is not necessary to have a R-parity violating contribution. 
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I. INTRODUCTION 



Neutron and electron electric dipole moments (EDMs) are important topics for new CP 
violating physics. They are known to be extremely small in the Standard Model (SM); in 
fact, way below the present experimental limit. With supersymmetry (SUSY) comes many 
plausible extra EDM contributions. That has led to the so-called SUSY CP problem |jj 
for the minimal supersymmetric standard model (MSSM). If one simply takes the minimal 
super symmetric spectrum of the SM and imposes nothing more than the gauge symmetries 
while still admitting soft SUSY breaking, the generic supersymmetric standard model would 
result. When the large number of baryon and/or lepton number violating terms in such a 
generic supersymmetric SM are removed by hand, through imposing an ad hoc discrete 
symmetry called R parity, one obtains the MSSM Lagrangian. In the case of R-parity 
violation, two recent papers focus on the contributions from the extra trilinear terms in the 
superpotential and conclude that there is no new EDM contribution at the 1-loop level [Q. 
Perhaps it has not been emphasized enough in the two papers that they are not studying 
the complete theory of SUSY without R parity, which is nothing other than the generic 
supersymmetric SM; in particular, they have nelgected admissible RPV parameters other 
than the trilinear ones in the superpotential. It is interesting to see that in the generic case 
there are in fact contributions at the 1-loop level, as pointed out in Refs. @,[§]. In particular, 
Ref. [0 gives a clear illustration of the much overlooked existence of a R-parity violating 
(RPV) contribution to LR squark mixings and the resulting contribution to neutron EDM 
through the simple 1-loop gluino diagram. We would like to emphasize again that the 
new contribution involves both bilinear and trilinear (RPV) couplings in the superpotential. 
Since RPV scenarios studied in the literature typically admit only one of the two types 
of couplings, the contribution has not been previously identified. A simple estimate of 
the bound obtained on the RPV parameters (the ji* X in combination) given in Ref. || has 
already illustrated that the bound from the neutron EDM as one of the most important, 
being competitive even when compared with sub-eV neutrino mass bounds and substantially 
more stringent than most collider bounds. The present article aims at giving a detailed 
analysis and numerical study of the RPV extension of SUSY contributions to neutron EDM. 
Similar new RPV contributions to electron EDM have been noted in Ref. 0. In fact, the 
complete result for RPV contributions to the masses of the sleptons and other (color-singlet) 
scalars has been given in Ref. [g], which focuses mainly on their implications to neutrino 
masses. 

The complete theory of SUSY without R parity admits all kinds of RPV terms without 
bias. It is obviously better motivated than ad hoc versions of RPV theories. The large 
number of new parameters involved, however, makes the theory difficult to analyze. The 
question of the specification of flavor bases to define the parameters in the Lagrangian of 
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the theory unambiguously becomes more important. In fact, thinking about the theory 
as the generic supersymmetric SM instead of as "MSSM + RPV terms" helps to clarify 
many of the issues involved ||. From such a perspective, it has been illustrated []/]] that an 
optimal parametrization, called the single- VEV parametrization (SVP), provides a very nice 
formulation which helps to simplify much of the analysis. In particular, the SVP gives the 
complete results for the tree-level mass matrices of all state, fermions as well as scalars, in 
the simplest form |J. The formulation has been used to study leptonic phenomenology J7] 
and various aspects of neutrino masses |5|JHTlll- The present EDM study (also Refs. [H,|J) 
and parallel works on /i -> e 7 (see also Ref. |13[|), electron EDM, and 6^57 ]I4|] will 



further illustrate the advantage of adopting the SVP. 

We focus here only on such contributions to the neutron EDM, based on the valence 
quark model ]T5|] . Hence, we study only the 1-loop quark EDMs. We will give complete 
1-loop formulae for EDMs of the up- and down-sector quarks, of which the u and d results 
are used to calculate the neutron EDM through the 

d n = - {Add- d u ) rj (1) 

formula, where rj ~ 1.53 is a QCD correction factor from renormalization group evolution 
rSyTTfl . This is to be matched with the experimental bound [|T3[ 



d n < 6.3 ■ i(T 26 e ■ cm . 

In the MSSM case, one has the SUSY loop contributions and the charged Higgs contribu- 
tions. The latter are very negligible. We focus here in this article on the analogue of the 
1-loop SUSY contributions. The latter include the gluino loop, the charginolike loop, and 
the neutralinolike loop. By the last two, we mean generalization of the chargino and neu- 
tralino loops under the generic picture. The (RPV) mixings of the leptons with the gauginos 
and Higgsinos give five (color-singlet) charged leptons and seven neutral fermions, including 
the charginos and neutralinos as well as e, fi, r, and three physical neutrinos. They come 
from the same set of electroweak states and should not be separated from one another in 
the analysis. It is no surprise that the physical chargino and neutralino states dominate 
the EDM contributions. We use explicit exact mass eigenstate expressions in our analysis 
to illustrate that as well as other interesting features, starting from the generic electroweak 
states couplings under the SVP. An exact numerical calculation is performed. We would 
like to mention that the generalization of the charged Higgs loop contribution involves other 
different RPV parameters. Moreover, there are many new and potentially important con- 
tributions including a t-quark loop, as also pointed out in Ref. [|J. We will give also the 
formulae of such quark-scalar loop contributions, though a detailed study is postponed to a 
future publication. Note that Ref. ||, which first appeared around the same time as Ref. 
||, is the only other study of the same topic available. In our opinion, our study here is 



3 



more systematic and complete. Ref. |4j] does not include, for example, the RPV LR scalar 
mixing and the resulted gaugino loop contribution to EDM. Moreover, to the best of our 
knowlegde, the present study includes the first exact numerical calculation performed. Ref. 
Q also quotes a cos/3 dependence of the major charginolike contribution, hence a weakening 
of the bound in the large tan/? regime — a result with which we disagree. Our careful nu- 
merical study illustrates many more interesting features, as the discussion below will speak 
for itself. 

This paper is organized as follows: We first summarize the formulation and notation used 
in Sec. II, where we also elaborate in some detail on the electroweak fermion field couplings 
needed to study the quark EDMs. Next, Sec. Ill contains results presented recently in Ref. 
H on all the scalar masses, listed here so that the present paper will be self-contained. Of 
most importance here is the RPV contributions to LR mixings, which play a central part 
in the EDM contributions. The quark EDMs are analyzed in Sec. IV. Some results from 
our numerical calculations are presented in Sec. V, after which we conclude in Sec. VI. An 
appendix gives some background formulae on the color-singlet fermion masses. Note that 
our formulae and calculations here naturally include the R-parity conserving MSSM part, 
though we will concentrate on the role of the RPV parameters and their unique contributions. 



Corresponding studies of the MSSM case can be found, for example, in Refs. []T6| , [T9| pi]] , to 
which we refer the readers for comparison. 



II. FORMULATION AND NOTATION 

We summarize our formulation and notation below. The most general renormalizable 
superpotential for the generic supersymmetric SM can be written then as 



W = e ab 



fi a H"L b a + y^ k Q1H b u U k + X ajk L a a Q b D k + - \ a p k L a a L b pE k + - X' ijk Uf 'D^ ' D k , (2) 
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where (a, b) are SU{2) indices, (i,j,k) are the usual family (flavor) indices, and (a, (3) are 
the extended flavor indices going from to 3. At the limit where A^, X i j k , A" j7c , and \ii all 
vanish, one recovers the expression for the R-parity preserving case (i.e., MSSM), with L 
identified as H^. Without R parity imposed, the latter is not a priori distinguishable from 
the Lj's. Note that A is antisymmetric in the first two indices, as required by the SU(2) 
product rules, as shown explicitly here with e 12 = — e 21 = 1. Similarly, A" is antisymmetric 
in the last two indices from SU(3) C - 

R parity is exactly an ad hoc symmetry put in to make L stand out from the other 
Li's as the candidate for Hd- It is defined in terms of baryon number, lepton number, and 
spin as, explicitly, 1Z = (— i) 3 - B + L + 2S '. The consequence is that the accidental symmetries 
of baryon number and lepton number in the SM are preserved, at the expense of making 
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particles and superparticles having a categorically different quantum number, R parity. The 
latter is actually not the most effective discrete symmetry to control superparticle mediated 



proton decay , but is most restrictive in terms of what is admitted in the Lagrangian, or 
the superpotential alone. 

A naive look at the scenario suggests that the large number of new couplings makes the 
task formidable. However, it becomes quite manageable with an optimal choice of flavor 
bases, the SVP ]7). In fact, doing phenomeno logical studies without specifying a choice 
of flavor bases is ambiguous. It is like doing SM quark physics with 18 complex Yukawa 
couplings instead of the 10 real physical parameters. For SUSY without R parity, the choice 
of an optimal parametrization mainly concerns the 4 L a flavors. Under the SVP [] , flavor 
bases are chosen such that 1/ among the L^s, only L , bears a VEV (i.e., (Li) = 0); 
2/ Vjk(= X ojk = -Xjok) = ^ diag{m 1; m 2 , m 3 }; 3/ yf k (= X 0jk ) = & diag{m d , m s , m b }; 4/ 
Uik = ^ KScm diag{m M , m c , m t }, where v = y/2 (L ) and v u = y/2 (H u ). The big advantage 
of here is that the (tree-level) mass matrices for all the fermions do not involve any of the 
trilinear RPV couplings, though the approach makes no assumption on any RPV coupling 
including even those from soft SUSY breaking, and all the parameters used are uniquely 
defined, with the exception of some removable phases. 

We are interested in scalar-fermion-fermion couplings similar to those of the charginos 
and neutralinos in the MSSM. The gaugino couplings are, of course, standard. Coming from 
the gauge interaction parts, they have nothing to do with R parity. The "Higgsino-like" 
part is, however, different. Without R parity and in a generic flavor basis of the four L Q 's, 
the Hd of the MSSM is hidden among the latter. The SVP, however, identifies L as the 
one having "Higgs" properties of H^, though it still maintains (RPV) couplings similar to 
those of the Lj's. We write the components of a L a fermion doublet as 1° and l~, and their 
scalar partners as 1° and l~. Apart from being better motivated theoretically, the common 
notation helps to trace the flavor structure. However, we will also use notation of the form 
h* and h* as alternative notation for /* and l£ in some places below. This is unambiguous 
under our formulation. We will also referred to h* and h* as the Higgs boson and Higgsino, 
respectively, while they are generally also included in the terms slepton and lepton. 

Note that in the left-handed lepton and slepton field notation introduced above, we have 
dropped the commonly used L subscript, for simplicity. For the components of the three 
right-handed leptonic superfields, we use If and If, with again the R subscript dropped. 
The notation for the quark and squark fields will be standard, with the L and R subscripts. 



1 Note that our notation here is a bit different from that in Ref. 0; we follow mostly the notation 
in Refs. || and || while improving and elaborating further whenever appropriate. We will clarify 
all notation used as our discussion goes along ||. 
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A normal quark state, such as d L , denotes a mass eigenstate, while a squark state the 
super symmetric partner of one. A quark or squark state with a / denotes one with the quark 
state being the SU (2) partner of a mass eigenstate. For instance, u' L3 is the up-type squark 
state from Q 3 which contains the exact left-handed b quark according to our parametrization 
of the Lagrangian. 

The up-sector Higgs boson is unaffected by R-parity violation. The scalar and fermion 
states of the doublet are denoted by h+, W and h£, ft, respectively. 

The identity of the charginos and neutralinos is unambiguous in the MSSM. Without R 
parity, they mix with the charged leptons and neutrinos. In fact, the true charged leptons 
and neutrinos are the light mass eigenstates of the full 5x5 and 7x7 mass matrices, respec- 
tively. The mass eigenstates deviate from the l^s and l°'s. Though the latter deviations are 
practically negligible in the limit of small /ij's, it still helps to distinguish the electroweak 
states from the mass eigenstates. Moreover, large R-parity violation, especially in terms of 
a large /x 3 , is not definitely ruled out 0. Here, in this paper, we make the explicit distinc- 
tion and reserve the terms chargino and neutralino for the heavy states beyond the physical 
charged leptons e, /i, and r and neutrinos. The m^'s introduced above are the Yukawa con- 
tributions to the physical charged lepton masses, hence approximately equal to the latter. 
The readers are referred to the appendix for more details about the fermion mass terms. 

We are now ready to spell out the couplings concerning the (color-singlet) charged and 
neutral fermions from the superpotential. We have 



c x = v clm K [ u c R . d h . + u c R . d L . ] + y d . I' [ d c R , u' L . + d c R , u' L . ] + X ijk I 



u 



v£. u L , + uiu L 



+ Vh 



+ u 



dl, dr +di .d L 



V 7° 

A ijk l i 



do, dr . ~\~ d~ dr 



r loif _ r 70 ;+ 

•'0 1 b l I ''O 1 



1° 1+ 7" 



+ 



+ Vh 
h.c. , 



7?7+r _ 7° 7+r 

1 l ''0 ''o 2 1 



(3) 



where 



Vu, 



g 2 m Ui 



V2M W sin/3 



92 mcL 



92 mi 



V2M W cos/3 



y/2M w cos/3 



(4) 



are the (diagonal) quark and charged lepton Yukawa couplings, and tan/3 = ^ ||23|| . Recall 
that X j k corresponds to the down-quark Yukawa coupling matrix, and \ jk corresponds to 
the charged lepton Yukawa coupling matrix, both of which are diagonal under the SVP; 
in addition, we have u' L , = V^^u^. being the SU{2) partner of the mass eigenstate d L ., 
and u' L its scalar partner. We also use below d[,, which is, explicitly, V^ M d Lj . There are 
some more scalar-fermion-fermion terms besides those given in C x . These extra terms are 
slepton-quark-quark terms. We will see that the latter are actually also involved in 1-loop 
EDM diagrams, though not the major focus of this paper. With the above explicit listed 
terms, however, it is straightforward to see what the extra terms are like. 
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In both of the above expressions for C x , there is a clear distinction between the MSSM 
terms and the RPV terms. The nice feature is a consequence of the SVP. The simple 
structure of the trilinear coupling contributions to the (i-quark and charged lepton masses, 
which is equivalent to that of the R-parity conserving limit, is what makes the analysis 
simple and easy to handle. We want to emphasize that the above expressions are exact 
tree-level results without hidden assumptions behind their validity. They are good even 
when there is large R-parity violation. 



III. SQUARK AND SLEPTON MASSES 



The soft SUSY breaking part of the Lagrangian, in terms of scalar parts of the superfield 
multiplets, can be written as follows: 



V soit = H a u L b a + e« b \A% Q?H b u U? + A%H a d Q\^ + A^H^Ef 



+ e a6 



13 
1 



+ h.c. 



A> ijk L iQj D k + 2 A %k L i L j E k 



+ - A) jk U°D°D c k + h.c. 



+ Q ] ml Q + U^rrf U + D^m 2 D D + L ] m 2 L L + &ih 2 E E + m 2 H jH u \ 



(5) 



where we have separated the R-parity conserving ones from the RPV ones (Hd = Lq) for 
the A terms. Note that L^m 2 L, unlike the other soft mass terms, is given by a 4 x 4 matrix. 
Explicitly, m 2 oo is rh 2 d of the MSSM case while ^ ofc 's give RPV mass mixings. The other 
notation is obvious. 

The SVP also simplifies much the otherwise extremely complicated expressions for the 
mass-squared matrices of the scalar sectors. First, we will look at the squarks sectors. 
The masses of up squarks obviously have no RPV contribution. The down-squark sector, 
however, has an interesting result. We have the mass-squared matrix as follows: 



M 2 = 



M 2 LL M 2 [ 
M 2 RL M 2 RR 



(6) 



where 



M 2 L =m 2 + mlm D + M 2 cos2/3 
M 2 RR = fnl+ m D m ] D + M 2 cos2/3 



~ + ^sin 2 6t 
-isin 2 ^ 



(7) 



and 



2 NT 
RLJ 



V2 



. V* a X ajk ) -j= 



= [A d - //* tan/3] 77^ + ^k^l 5A ° 



\/2M w sin/3 



92 



92 



(^iKjk ) • 



(8) 
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j k with elements listed [] . The (A**A' ajA ) term is the 
{ j k j part separated out in the last expression gives 



Here, m D is the down-quark mass matrix, which is diagonal under the parametrization 
adopted; Ad is a constant (mass) parameter representing the "proportional" part of the A 
term and the matrix 5A D is the "proportionality" violating part; (/^A'^), and similarly 
( £**A' - fc )j denotes the 3x3 matrix 
full F-term contribution, while the 
the RPV contributions. 

Next we go on to the slepton sector. From Eq.(|5]) above, we can see that the charged 
Higgs bosons should be considered on the same footing together with the sleptons. We have 
hence an 8 x 8 mass-squared matrix of the following 1 + 4 + 3 form: 

( Ml, M& Mil \ 



M; 



M* Ml Ml\ 



M 2 

RH J RL 



(9) 



where 



M: 



x'2 



m 



4+/i> Q + M^cos2/? 



- sin 2 ft- 



+ Ml sin 2 /? [1 - sin 2 6U] , 
M 2 LL = ml + m L m L + + M| cos2/3 



— - + sin 2 8 w 



+ 



Ml cos 2 /? [1 - sin 2 6W 



1X3 




M' 



m E + m E m E + M z cos2/3 — sin 9 W 



3x3 
2/ 



(10) 



and 



ML 



M z 



(M 



2 NT 

rl) 




\ Ml sin2/? [1 - sin 2 ^ 



v ° i * \ 
-j= = {fx k m k ) 



(no sum over k) 



- (KA Q(3fc ) 



[A e - /x* tan/3] 




V2 

\/2 M w cos/3 




— fil m k tan/3 

y^ilfr sin/3 
92 



The notation and results here are similar to the squark case above, with some difference. 
We have A e and 5A E , or the extended matrices ^ ° ^ incorporating them, denoting the split- 



ting of the A term, with proportionality defined with respect to m, 



diag{0,m B } 



2 Note that we use this kind of bracket notation for matrices extensively here. In this case, the 
repeated index i is to be summed over as usual and, hence, is dummy. 



S 



diag{0, n\, rrh, m 3 }. Recall that the m/s are approximately the charged lepton masses. A 
4x3 matrix (fj,*X i/3k ) gives the RPV contributions to (MlY '. In the above expression, we 
separate explicitly the first row of the former, which corresponds to mass-squared terms of 
the type l + h~ type (h~ = The nonzero Ml and the _B*'s in Ml are also RPV con- 
tributions. The former is a l + (h^y type, while the latter a l~h+ term. Note that the parts 
with the [1 — sin 2 6^] factor are singled out as they are extra contributions to the masses of 
the charged Higgs bosons (i.e., Z~ = h~ and The latter are the result of quartic terms 
in the scalar potential and the fact that the Higgs doublets bear VEVs. 

The neutral scalar mass terms, in terms of the (1 + 4) complex scalar fields ra 's, can 
be written in two parts — a simple (-M 2 0t ) mri 0j„0 n part and a Majorana-like part in the 
form | (M^^mn (pm<Pn + h.c. As the neutral scalars originate from chiral doublet superfields, 
the existence of the Majorana-like part is a direct consequence of the electroweak symmetry 
breaking VEVs, hence restricted to the scalars playing the Higgs boson role only. They come 
from the quartic terms of the Higgs fields in the scalar potential. We have, explicitly, 



2 z 



( sin 2 /3 — cos/3 sin/3 



— cos/5 sin/5 



cos 2 /3 



o lx3 \ 

0lx3 

0;ix3 / 



(12) 



and 



M = M + 



rhk + + M z 2 cos2/3 

-m 



-{Bo) 

rhl + to/0 + M z cos2/? 



(13) 



Note that M.^ here is real (see the Appendix), while M^ t does have complex entries. The 
full 10 x 10 (real and symmetric) mass-squared matrix for the real scalars is then given by 



M' 



Ml 



(Mir m. 




(14) 



where the scalar, pseudoscalar, and mixing parts are 



Ml = Re(M 2 J + M, 



Mi 



MM*) - Mi 



Ml = -lm(ML) , 



(15) 



respectively. If lm(Ml) vanishes, the scalars and pseudo-scalars decouple from one another 
and the unphysical Goldstone mode would be found among the latter. Finally, we note that 
the B a entries may also be considered as a kind of LR mixing. 

We would like to emphasize that the above scalar mass results are complete — all RPV 
contributions, SUSY breaking or otherwise, are included. The simplicity of the result is 
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a consequence of the SVP. Explicitly, there are no RPV A-term contributions due to the 
vanishing of VEVs Vi = y2(Li). The Higgs-boson-slepton results given as in Eqs.(§) and 
(|TJ) contain a redundancy of parameters and hide the unphysical Goldstone state. However, 
the results as they are given here are good enough for some purposes including the present 
EDM discussion. 

Note that in the following discussion, we will not consider flavor changing scalar mass 
terms from soft SUSY breaking, which could be suppressed through a flavor-blind SUSY 
breaking mediating mechanism such as gauge mediation. A major concentration of our 
study, however, will be on the effects of the flavor changing scalar mass terms from RPV 
superpotential parameters, which give interesting new results. 

IV. CONTRIBUTIONS TO NEUTRON EDM AT 1-LOOP 

In this section, we will discuss the RPV 1-loop contributions to neutron EDM systemat- 
ically, drawing comparison with the R-parity conserving (MSSM) case wherever it would be 
useful. We will give complete 1-loop formulae for quark EDM's for our generic supersym- 
metric SM. We will not, however, go into the chromoelectric dipole operator or Weinberg 
gluonic operator contributions. Following the common practice, family mixing will largely 
be neglected, though we will comment on some particularly interesting aspects of the issue. 
We will also compare the results with Feynman diagrams given more or less in an electroweak 
basis. Naively, such diagrams, with a minimal number of mass insertions admitted, should 
represent a first order result, at least where mass mixings are small. Note that the latter is 
true for the small-/ij scenario, which is our major focus. Afterall, a mass eigenstate result 
may be considered as corresponding to taking an infinite summation over electroweak-state 
diagrams with all possible mass insertions. The comparison helps to illustrate better the 
physics hidden behind the complicated formulae. In addition to a short Letter by the present 
authors [[|, some parts of the results here have also been discussed in Ref. [|J. 

A. Gluino Contributions 

The most direct contributions come from a gaugino loop, as shown in Fig. 1. The diagram 
looks the same as the MSSM gluino and neutralino diagrams with two gauge coupling 
vertices. As pointed out in our previous short Letter ||, the new RPV contributions here 
are a simple result of the RPV LR squark mixings [cf. Eq.(|)]. In Ref. H, we focused 
on the dominant gluino loop contribution. We first give some more details of that analysis 
before going into the other contributions. Notice that though both the u and d quarks get 
EDMs from gaugino loops in the the MSSM, only the d quark has the RPV contribution. 
The w-squark sector simply has no RPV LR mixings. In Fig. 1, as well as the subsequent 
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diagrams, we use a family index k for the external quark lines though only the k = 1 case 
would be the d- (or u-) quark EDM we are mainly concerned with. With the only possible 
exception of d 2 , which corresponds to the s quark [[T5[], k 7^ 1 results are not relevant for 
neutron EDM. 

For illustrative purposes, we first neglect interfamily mixings among the squarks. The d 
mass-squared matrix, of the form given in Eq.(^) but reduced to one family, is Hermitian 
and can be diagonalized by the unitary transformation 



VlM 2 D V d = dmg{Ml,Mi + }, 



with 



and 



Ml 



(M 2 LL + M 2 RR ) =f J (Ml - Ml R Y + 4 \M h 



Here, M\ h = 
tanfl - 



-Mi 



(16) 
(17) 

(18) 

|-M|J e 1 ^ and the range of 9 can be choosen so that — | < 9 < | with 
■. In terms of the mass eigenstates gL and d + , the gluino contribution to 



cos ~ 
sin | e -< * 



sin | e -< * 
cos ~ 



the EDM is then given by JT|,|20 

im(r d 





2« s 


(tI- 


3tt 






where Q d ~ = — |, 


rife _ 



a— 



Mr 
a- 



+ im(r, 



12 1 M§ O HI 



M 2 



(19) 



V d2k V* dlk , giving Im(r d 
1 



ir 



fix 



2(x- l) 2 



1 + x + 



-lm(Tj 2 ) = i sin9 sin0, and 
2x\nx 



x 



Introducing M| = (M£ + M? + )/2 and using the identity relation xB(x) 
gluino contribution becomes the often-quoted 



— 2 QjIm^jF 



,/M 2N 



I e /- 3tt M 2 

\ / g d 

where 8° is M. RL /M? (with .M 2 ^ restricted to the g? family) and 

Fix) 



(20) 
the 

(21) 



'i-xY 



1 + 5x (2 + x)x In a; 



1-x 



(22) 



The EDM expression above is, in fact, the same as that of the MSSM case, except that 
5° , or equivalently A4 RL , has now an extra RPV part. From the general result given in 
Eq.(||), we have, for the d squark, 
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[A d - ji* tan/3] m d 



V2M w cosB AD V2M w sm(3 

~ ° A u : ( IH \n 



9> 



(23) 



Note that the /i*A' in term does contain nontrivial CP violating phases and gives RPV con- 
tributions to neutron EDM. Though the above analysis neglects interfamily mixings, it is 
clear from Fig. 1 that the LR squark mixing 8° is what is directly responsible for the EDM 
contribution. Including interfamily mixings would complicate the mass eigenstate analysis 
but not modify the EDM result in any substantial way. 

B. Neutralinolike Contributions 

The contributions from the (electroweak) neutral gaugino loop, as shown in Fig. 1, are 
expected to be small, due to the much smaller gauge couplings. Apart from the neutral 
gaugino loops, there may be other neutralinolike loop contributions. In the MSSM, one 
has to consider possible contributions from the Higgsino parts and the gaugino- Higgsino 
mixing parts. Part of such contributions involves no LR squark mixing. The latter feature 
compensates for the smaller Yukawa-type couplings involved [|19j. Without R parity, the 
gauginos and Higgsinos mix with the leptons. We have seven neutral fermions, all among 
which give rise to a neutralinolike loop contribution. 

The neutral fermion loop contribution to u and d quark EDMs is given by 
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with b = 3 (4) for T 3 f = ~ (—5) [i.e., for / being the u (d) quark], and y f the corresponding 
Yukawa coupling [cf. Eq.(f|)]. Recall, from Eq. (|16|) , that d T denotes the two d mass eigen- 
states and T> d the diagonalizing matrix; likewise, and T> u are the corresponding notations 
for the w-quark case. Finally, the -X^-'s are elements of the matrix that diagonalizes 
as defined explicitly in the Appendix. Note that the last term in each set of brackets of the 
A// n=F expressions is nonvanishing only for / = d, as indicated by the 8^ symbol. Similar 
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to the case of the gaugino contributions, RPV contributions here exist only for the <i-quark 
EDM. We would like to emphasize that the formulae here represent the full 1-loop result for 
our generic supersymmetric SM. 

Each expression for A// nT above is a product of two terms, from the two loop vertices 
involving the L- and -R-handed quark fields, respectively, as a chirality flip has to be imposed 
within the loop. For each vertex, the three terms (in each set of brackets) correspond to 
the gauge, quark Yukawa, and RPV A' couplings, with the last existing only for the (i-quark 
case. With two gauge coupling vertices, we have the gaugino diagram. Diagrams with a 
gauge and a y f coupling are shown explicitly in Fig. 2 for the d quark, in which y d is indicated 
by X okk , i.e., with the notation used in our superpotential and a general, unspecified, quark 
family index ( y d = X on in our notation ). Note that the diagrams require no LR mixing 
on the squark line. This is the important MSSM contribution, the chargino counterpart of 
which typically receives the major attention. The latter is numerically a bit larger. The 
type of diagrams has no RPV contribution. 

There is, however, a RPV analogue to Fig. 2. With the notation as given, this is obviously 
given by replacing X okk in the figure with a RPV X ikk . This is shown explicitly in Fig. 3. 
From the latter figure, one can see the the first order result would come from a /"-gaugino 
mass mixing. However, under our formulation, the latter is vanishing [see Eq. (|A8|) in the 
Appendix]. Looking at the type of contributions from our EDM formula here, while a 
X kn coupling may not be small, it comes into the formula with a X^+4)n- For n = 1-4, 
corresponding to the heavy neutralino mass eigenstates, X^ k+4 ^ n measures a RPV mixing 
in the neutral fermion masses M.^- For n = 5-7, at least one of the three X(k+4)nS is 
expected to be of order 1, but the corresponding (physical neutrino) mass eigenvalues M^o's 
give a strong suppression factor in the resulting EDM contributions. There is also a further 
suppression from the mixing factor of the gaugino part, e.g., Xi(k+4)- Furthermore, there is, 
potentially, a GIM cancellation among the seven mass eigenstates. Our numerical results, 
however, do show that the type of contribution is generally important. We will discuss the 
issue more carefully below, using the explicit example of its charginolike counterpart. 

Last, we come to the diagrams with no gauge vertices, which we show in Fig. 4 for the d 
quark, using our generic A' fe notation. Such a diagram in the MSSM case is totally negligible 
due to the small Yukawa couplings involved and the suppression from the LR squark mixing 
required. In fact, as shown in the figure, the minimal mass insertion needed for the diagram 
corresponds to a Majorana mass among the /°'s, which is again vanishing. The situation is 
similar to that of Fig. 3. The contribution is expected to be always dominated by the latter 
one. 
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C. Charginolike Contributions 



Next we come to the charged fermion counterpart. With similar notation as used in 
the neutral fermion loop formula above, the charged fermion loop contributions to u- and 
d-quark EDMs can be written as 
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for / being u (d) quark and /' being d (u), where 
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(only repeated index i is to be summed) ; (27) 

V and U are diagonalizing matrices of the charged fermion mass matrix M. c as defined by 
Eq.(^) in the Appendix; and the function A(x) is given by 

w s 1 21nx\ 

The first and second parts of Eq.(^) come from the cases in which the photon is emitted 
from the squark and fermion lines inside the loop, respectively. 

The basic feature of the C/ n=F 's terms is similar to the previous case of A// n=F . They do 
not give rise to a charged gaugino loop contribution, though, as i?-handed quarks do not 
couple to W ± . Taking the available gauge coupling term within each C/ n=F and a y f to form 
an admissible product, we have the contribution corresponding to a diagram in Fig. 5. This 
is the dominating MSSM contribution, besides the equally competitive gluino one. Note 
again that no LR squark mixing is required the diagrams. 

The diagram given in Fig. 5b has the RPV analogue, which requires a l~ k -W + mass 
insertion for the first order result, as shown in in Fig. 6. This is the SU(2) partner of Fig. 3, 
something we promise to discuss in more detail. Looking at the contributions from our EDM 
formula here, we have again a X kn coming in with a U^ +2 ) n , as versus the X( k+4 ^ n in the 
neutral case above. For n = 1-2, corresponding to the heavy chargino mass eigenstates, 
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U(k+2)n measures a RPV mass mixing; it is of order 1 for n = k + 2, but then M- is the 
small nik, roughly the charged lepton mass. When one sums over for n = 1-5, it is easy to 
see that the result is proportional to the imaginary part of 

£ V* ln M x -U (k+2)n F SA (M 2 .) A' feu , (29) 

71=1 

where F BA iM 2 _ J denotes a function on M\ corresponding to the large brackets in Eq.(p6|) 
with functions B and A. If the F BA could be factored out, together with A' fcu , what is left inside 
the summation is nothing other than the l k -W + mass term, which is zero. This is a GIM-like 
cancellation mechanism, violated only to the extent that the loop integrals involved as given 
by the B and A functions are not universal. Our numerical calculations, however, show 
that the cancellation is very substantially violated, for generic values of chargino masses. 
Let us then look at the contribution from an individual mass eigenstate more closely. To 
get an analytical approximation, we used the perturbative results for U^+2) n given in the 
Appendix. For the n — 1 and 2 parts in the above sum, we have 

V* la ^ k R R2a F BA (Ml) X kn , (30) 

with basically the same source of RPV complex phase as the gluino case, namely, from 
lm(filX kn ), except that it is one value of k for each diagram here. Note that the explicit 
M Ca factors are canceled. The n = k + 2 part is largely suppressed due to a small mass 
eigenvalue and the very small RPV i?-handed mixing given by VJV fc+2 s as shown in the 
Appendix. Note that Vj* a ~ R£ , and R R is a 2 x 2 unitary matrix of order 1 elements. 
The expressions give an idea about the strength of the RPV contributions. In the limit of 
degenerate charginos, i.e., M cl = M c2 , the F BA function factors out of the sum of the n = 1 
and 2 parts and the GIM-like cancellation is clearly illustrated, simply from the unitarity 
of R R . In fact, our numerical calculations give a cancellation up to 1 part in 10 4 in such a 
situation if we enforce the condition. However, that requires a every substantial complex 
phase for /x , hence of not the most phenomeno logical interest. Another interesting point on 
the parameter space appears when the R R matrix is diagonal, at which the contribution also 
goes away, as is obvious from the above expression. This happens at the point where the 
condition M* sin/3 = fi cos/5 is satisfied. This feature will be confirmed by our numerical 
results presented below. 

The rest of the charginolike contributions each has at least one y f coupling and no 
gauge coupling vertex. Again, a LR squark mixing is required. However, of the first order 
electroweak-state diagrams, as depicted in Fig. 7, each has an admissible /z* mass insertion 
on the fermion line. This looks a bit different from the neutral case above. Apart from 
an admissible RPV down-squark mixing for the case of the u quark, a RPV diagram again 
has the RPV parameter combination ^ k X kn (no sum). Except for large there is no 
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reason to expect the contribution to be of any significance compared to the other major 
contributions. 

Our EDM formulae neglect interfamily mixings. However, general, unspecified, family 
indices are use in our figures to make the flavor structure transparent. We note that in the 
case of the contributions depicted by Figs. 4 and 7, off-diagonal /^A'^-type RPV mixings 
could play a role to let the higher family squarks mediating more important contributions 
to u- and d-quark EDMs. 



D. The Quark-Scalar Loop Contributions in Brief 



There are superpartners to the diagrams in Fig. 7, i.e., with quarks and charged scalars 
running inside the loop and a B a coupling on the scalar line. These diagrams are shown 
explicitly in Fig. 8. Here, the a = case gives the MSSM charged Higgs boson contribution. 
With the -Bj's, we have the RPV contributions. Analytically, we have the formula 
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and, for f = d, 



nR _ 



92 

Ct m = ^V l 2m + ^V\ t+2)m , (33) 
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with the J2' m denoting a sum over (seven) nonzero mass eigenstates of the charged scalar; 
i.e., the unphysical Goldstone mode is dropped from the sum, V 1 being the diagonalization 
matrix, i.e., V^M 2 E V l = diag{ A^,m = 1-8}. 

With similar notation as used in the charged scalar loop formula above, the neutral scalar 
loop contributions to the EDMs can be written as 




, 47rsin^ tT^i V mm mmJ Ml ^ \Ml 

where, for / = u, 
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with again the unphysical Goldstone mode to be dropped from the sum over scalar mass 
eigenstates (nine nonzero), and the obvious notation {T> S ) T M^ V s = diag{ M? , m = 1 - 10 }, 
V s being an orthogonal matrix. 

Note that the formulae given in this section, like those in the rest of the paper, have 
neglected CKM mixings among the quarks. However, unlike the previous cases, the mixings 
play an important role in the quark-scalar loop contributions. This is a result of the fact 
that the EDM contributions have a fermion, a quark in this case, mass dependence, and 
the mass hierarchy among the quarks. For instance, the top quark loop is expected to give 
a dominating contribution. Generalizing the formulae to include the mixings is straightfor- 
ward. The formulae as they are given in this subsection do give a basic illustration of the 
major analytical structure of the type of contributions. We will only discuss the features 
briefly here. 

Diagrams depicting the neutral scalar contributions are like superpartners of the type of 
diagram given in Fig. 3. This is shown in Fig. 8 for both the u- and <i-quark cases. Such a 
diagram could be interpreted as requiring a Majorana-like scalar mass insertion. The latter 
is first considered in Ref. ||11|| , under the name of Majorana masses from the "sneutrinos" . 
However, looking at it from the present framework, these diagrams should be interpreted as 
having Majorana-like mass insertions among the Z°'s (and W u for the case of the u quark). 

The quark-scalar loops given by the above formulae are not quite considered in the 
case of the MSSM, as they would be suppressed by the very small Yukawa couplings. The 
RPV analogue, however, admits much more general flavor structure, as illustrated in the 
diagrams in Figs. 8 and 9. For instance, a top quark with order-1 Yukawa coupling could be 
contributing to the <i-quark EDM through a charged scalar loop. Contributions discussed 
in this subsection depend on a much larger number of parameters, through the scalar mass 
matrices [cf. Eqs.([|)-(|r|)]. The extra parameters are related directly to Higgs physics. The 
important RPV Bi parameters here have a strong connection with the /Vs ||23|| . Furthermore, 
the more general flavor structure involved means that the EDM contributions involve 
couplings that would also have important roles to play in the related processes of b — ► s 7 
and b — > d 7. We will leave all these issues for later studies. 
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V. RESULT FROM NUMERICAL CALCULATIONS 



We now discuss the results we obtained by a careful numerical implementation of our 
EDM formulae discussed above, with explicit numerical diagonalization of all the mass 
matrices involved. Note that the quark-scalar loop contribution is not included in the 
numerical study. Hence, we are concentrating on the SUSY contributions in the presence of 
RPV couplings, to be compared directly with the MSSM results. 

As discussed in the previous section, the imaginary part of the combination fi*X in is what 
is central to the RPV 1-loop EDM contributions. To simplify the discussion, we single out 
fM 3 and A' 3U and put the corresponding parameters for i — 1 and 2 to be essentially zero []. 
This applies to all the numerical results discussed here. 

We give in Table 1 some detailed numerical results for four illustrative sample cases. To 
first focus on the RPV contributions which we are particularly interested in, cases A and B in 
the table have all complex phases in the R-parity conserving part switched off. For the RPV 
parameters /i 3 and A' 311 , we take the former being real and put a 7r/4 complex phase into the 
latter. However, we want to emphasize again that only the phase of the combination /x*A 31L 
matters. For instance, we have checked explicitly that shifting the overall phase, or a part 
of it into /i 3 instead, gives exactly the same results. The difference bewteen case A and case 
B is only in the value of /x 3 chosen. Case A has small /x 3 , at 10~ 3 GeV. This is the small-/^ 
scenario, corresponding to a sub-eV mass for v T as suggested, but far from mandated, by 
the result from the Super-Kamiokande (super-K) experiment [|23J []. In contrast, case B has 
fi 3 = 1 GeV. Note that imposing the 18.2 MeV experimental bound [25|] for the mass of v T 
still admits a relatively large /z 3 , especially for a large tan/3. Reading from the results in 
Ref. @, the bound is ~ 7 GeV at tan/3 = 2 and ~ 300 GeV at tan/3 = 45. As for X 3n , the 
best bound on the (from r — > nu) is around 0.05-0.1 p6fl . Hence, case B is still within 
the admissible region of RPV parameter space, though beyond the more popular small-/ij 
scenario. Moreover, we have illustrated in Ref. that the A 3U bound is not strengthened 
even when the above-mentioned stringently interpreted neutrino mass bound from super-K 
is imposed on the 1-loop neutrino mass contribution obtainable from the parameter. We 
therefore use the same A 3U magnitude of 0.05 for both, and actually all, cases in the table. 

Apart from the overall EDM results, we are interested in understanding the relative 



3 In the actual calculation, tiny but nonzero values are used to avoid possible problems of numerical 
manipulations. 

4 Note that /j, 3 cos/3 ~ 10 -4 GeV is actually expected, allowing a larger // 3 in the case of large tan/3 
(see, for example, Ref. M). Our choice of fj, 3 value here is, admittedly, quite pushing the limit. It 
is mainly for illustrative purposes. 
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values of the different pieces of the contributions discussed in the previous section. Hence, 
we show, in the table, the contribution from individual pieces, identified by the couplings of 
the two loop-vertices given inside the first column. They can be matched easily with terms in 
the formulae. The corresponding Feynman diagrams are marked whenever they are available 
inside the next column. Column 3 of the table indicates whether a LR squark mixing is 
involved in the contribution, and if the RPV part of the latter is available and needed for a 
RPV contribution to come in. The individual contribution that is vanishing under cases A 
and B corresponds to a diagram which does not admit a RPV contribution. This is the case 
for most entries for the w-quark EDM and the entries for the charginolike loop contributions 
to the (i-quark EDM not involving a A' coupling loop vertex. Note that a charginolike loop 
for the d quark involves LR mixing of u which has no RPV part. As emphasized in Ref. |§, 
the RPV contributions are far more prominent for the d quark part. Again, A' 3U is the only 
nonzero A' coupling introduced, though we use a more general A' notation within the table 
to remind the readers of the more general flavor structure admissible, as illustrated in the 
previous section. We also recall from the discussions there that a A' ijfe coupling always comes 
into the EDM picture accompanied by a RPV fermion mixing matrix element proportional 
to //*. The EDM results are basically in direct proportion to \i\ X in over the whole region 
of parameter space studied. Hence, comparing cases A and B, we see an exact factor of 
10~ 3 suppression in case A for contributions with one A' loop vertex or involving a RPV LR 
mixing, except the very small term. 

The overall EDM numbers of case A are still below the present experimental bound. This 
is smaller than a naive estimated result, as given for the gluino contribution by Eq. (f2T|) , for 
example, due to some unavoidable partial cancellation. For the gluino result, in particular, 
the cancellation is between the contributions from the two squark mass eigenstates, which 
would in fact be exact when the latter states are degenerate. This slightly weakens the EDM 
bound of the previously estimated Im(/z*A' ilL ) < 10 _6 GeV result given in Ref. ||. However, 
we see here that with a unification-type relationship imposed on the gaugino masses, we are 
likely to have a chargino contribution larger than the gluino one. Moreover, the values of the 
other SUSY parameters chosen here can be pushed to increase the EDM number, as discussed 
below. As for case B, the numbers are beyond the experimental bound, indicating that the 
RPV parameters, or the involved complex phase, would have to be further constrained. 

Case C of Table 1 gives a different scenario. Here, we allow complex phases, and hence 
EDM contributions, of the MSSM part. The RPV parameters are set real. This ensures no 
RPV EDM contribution from diagrams involving LR squark mixings. The interesting point 
to note here is that there is indeed nonzero RPV contributions from diagrams with RPV 
loop vertices. In particular, the RPV chargino contribution to <i-quark EDM is comparable 
in magnitude to its R-parity conserving counterpart. Hence, the existence of nonzero RPV 
parameters, even real, would change the EDM story of the MSSM. In other words, in the 
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presence of a complex /x , the neutron EDM bound actually constrains the magnitude of the 
the combination of RPV parameters given by /i*A' in , instead of just the phase of it. This is 
an important feature that has not been pointed out before. The result of the case serves, 
otherwise, as a check against other neutron EDM studies of the MSSM. 

Finally, we illustrate in case D a specific situation with both R-parity conserving and 
violating phases. Here, we pick a case of large tan/3, for which the larger /x 3 is still within 
the limit admitted by the stringent interpretation of the super-K bound ||. The negative 
phase chosen for X 3n gives a cancellation between the two contributions to the imaginary 
part of LR d squark mixing. This directly resulted in suppression of the EDM contributions 
involving the latter, such as the <i-quark gluino loop. The very interesting point to note here 
is that the same cancellation results in the chargino and neutralino loop Q , as explicitly 
illustrated by the entries of the terms with loop vertex couplings given by g ■ y d and g ■ 
X in in both cases. Essentially, when we have a small imaginary part for n* a X an (recall 
that a = to 3, A' on = y d ), e.g., an internal cancellation among the summed terms, the 
corresponding d quark EDM contributions, involving the suppressed LR mixing [ cf. first 
line of Eq.(H) ] or otherwise, are all suppressed. The chargino and neutralino contributions 
actually reflect some proportionality to the (F-term) LR mixings in a way similar to the 
gluino contribution. Note, however, that the A-term phase contributes to the gluino diagram 
while leaving the charginolike and neutralinolike diagrams not quite affected, and hence 
spoils the simultaneous cancellation achieved in the case D results here. Besides, as the 
w-quark EDM results do not have a complete matching analogue between the R-parity 
conserving and violating contributions, the simultaneous cancellation mechanism does not 
exist there. In the current case, the u-quark parts are suppressed due to large tan/3. 

The above sample cases illustrate some of the interesting features about the RPV con- 
tributions to neutron EDM. Next, we discuss how the EDM result varies with some basic 
parameters. Fig. 10 shows a logarithmic plot of (the magnitude of) the RPV neutron EDM 
result for value between ±2000 GeV, with the other parameters set at the same values as 
case A in Table 1. Note the the small (< 70 GeV) region has been ruled out for having 
too light a chargino mass As noted in the previous section, the gluino and chargino 
contributions compete with one another with one of them being dominating in a region of 
the parameter space. The dip on the C line for the charginolike contribution corresponds 



5 A chargino (neutralino) contribution means exactly the one with a physical chargino (neutralino) 
running in the loop. Hence, it is only a part of the fermion mass eigenstate sum in the analytical 
formulae given for the (total) charginolike and neutralinolike loop contributions. The latter, as 
illustrated in the previus section, are generally dominated by the chargino and neutralino mass 
eigenstates. 
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to the case of vanishing i?-handed mixing among the charginos, i.e., when the condition 
M 2 sin/5 + /i cos/? = is satisfied, as noted above, where the contribution essentially van- 
ishes. Note that in the region to the left of the point, the contribution becomes negative. 
The G line is horizontal, as the gluino diagram result has no dependence on /Jq. 

Next, we check the tan/? dependence and compare with the MSSM result. In Fig. 11, 
we give again a plot of the (total) EDM result for the same set of input parameters as 
in case A of Table 1, while varying tan/3. This is the line marked as "RPV only". The 
numerical result confirms our earlier discussion that there is not much sensitivity in tan/? 
[cf. expression pOj)], as indicated by the flatness of the line. This is in good contrast to 
the MSSM result, also based on the same set of input parameters except with nonvanishing 
phases for A and /v The third line, marked by "GSSM", gives the total result obtained 
from our formulae given above with the same nonvanishing phases for A and fa, as well as 
A' 3U . Note that the GSSM result here is more than the sum of the two other lines, due to 
the presence of RPV contribution even in the limit of a real A' 3U , as discussed for Case C of 
Table 1 above. The tan/? dependence, or the lack of, illustrated in the figure is quite generic, 
in a wide region of the parameter space. 

The dip in the GSSM line of Fig. 11 represents a point where the overall contribution is 
small as a result of the cancellation among the different pieces, a feature that has attracted a 
lot of attention lately in the case of the MSSM (see, for example, Ref. and Refs. [pO| ,pl[] ) . 
The MSSM line in the figure has a dip only beyond the range of the tan/? value shown. As 
also pointed out in a little bit different setting (with vanishing A phase) for Case D of Table 1 
above, the new RPV contribution modifies the overall picture and provides a possible new 
cancellation mechanism. The cancellation feature is better illustrated in Fig. 12, in which 
the variation against the A-term phase is shown, for the small and large tan/? cases. Again 
our GSSM result is compared with that of the MSSM. One can see that the presence of the 
RPV contribution shifts the cancellation points substantially. Finally, we show variations 
of the result as a function of the gaugino mass parameters here represented by M 2 . This is 
given in Fig. 13 with the four lines marked and correspond to the four cases of Fig. 12, each 
with the A-term phase set at the position of the dip as given in the latter figure. 



VI. CONCLUSION 

We have given explicit formulae and detailed discussions on the full 1-loop contributions 
to quark EDMs for the generic supersymmetric SM (without R parity). The extra, RPV, 
contributions are interesting additions to the R-parity conserving part. We have given 
results from an exact numerical study, illustrating various novel aspects. Our formulation 
emphasizes the universal structure of the R-parity conserving and violating parts, which also 
is illustrated itself well in the results. The experimental bound on the neutron EDM is hence 
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established as an important source of constraints on the model parameter space, including 
the RPV part. The 1-loop RPV contributions always involve the particular combination of 
RPV parameters given by n*X in , with little sensitivity to the value of tan/5. So far as the 
RPV parameters are concerned, this combination is well constrained by the EDM bound. 
This applies not only to the complex phase, or imaginary part of, the combination. Real 
fJ.*X in contribute in the presence of complex phases in the chargino and neutralino mass 
entries. Studies with either fj, r or A'- type couplings assumed to be zero miss the class of 
very interesting phenomenological features of SUSY without R parity. 
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APPENDIX A: NOTES ON THE FERMION MASSES 



Under the SVP, the (color- singlet) charged fermion mass matrix is given by the simple 
form 
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with explicit bases for right-handed and left-handed states given by (— iW + , ft, t[,Q, 1% ) 
and (— iW~, l~ , l~, l~ , ll ), respectively Here, we allow M 2 and all four \i a parameters to 
be complex, though the m^'s are restricted to be real, for reasons that will become clear 
below. Obviously, each fii parameter here characterizes directly the deviation of the /" from 
the corresponding physical charged lepton = e, fi, and r), i.e., light mass eigenstates. 
For any set of other parameter inputs, the mj's can then be determined, through a simple 
numerical procedure, to guarantee that the correct mass eigenvalues of m e , m M , and m T are 
obtained — an issue first addressed and solved in Ref. 0. The latter issue is especially 
important when /Vs not substantially smaller than /z are considered. Such an odd scenario 
is not definitely ruled out ]7). However, for the more popular small-/ij scenario, we have 
l~ £~, and deviations on m,j's from the (real) li masses are negligible. Note that the 
deviation of / 4 + from £^ is quite negligible in any case. 

We introduce unitary matrices V and U diagonalizing the R- and L-handed states with 



V ] M C U = diag{M- } = diag{M cl , M c2 , m e , m M , m T } . 



(A2) 



Here, the mass eigenvalues M- with n = 1 and 2, i.e., M cl and M c2 , are the chargino masses. 
Consider further 



V2M w sxap Ho 



(A3) 



with M° x and M° 2 being the chargino masses in the fa — limit. One can then write the 
diagonalizing matrices in the block form 



V 



R R -R R V^ 
V l va 



u 



R L -R L tf 

U L vn 



(A4) 



For fa M" a (a = 1 and 2), a block perturbative diagonalization can be performed directly 
on the matrix M e to obtain the following simple result 



23 



Uii+2)! ~ 




-Rfl2l 


U {i+2 )2 ^ 








rrti 




V(i+2)a - 


M c 





(a = 1 and 2) . (A5) 

Elements in the R R and i? L matrices are all expected to be of order 1. The above expressions 
illustrate the magnitudes of given matrix elements involving RPV mixings, as well as their 
dependence on the RPV parameters. We note that the L-handed mixings are roughly 
measured by the ratio of a Hi to the chargino mass scale, while i?-handed mixings are 
further suppressed by a charged lepton to chargino mass ratio, hence quite negligible under 
most consideration. 

Note that the notation here is different from that given in Ref. and many others 
in the literature. More explicitly, we have R- and L-handed mass eigenstates given by 
( x+n ) = V T [-iW + , h+, Z+, Z+ Z+ ] T and (*-„) = U ] [-iW~, l~, Z", Z", l; ] T , which form the five 
Dirac fermions 



Xr. 



A+n 



The C/(j +2 )a elements as given above show no obvious dependence on tan/3, though some 
nontrivial dependence is expected through the R R2a elements. Our exact numerical result 
also indicates a weak sensitivity on the tan/3 value here. On the other hand, we have, from 
Eqs.QAD and Q£5|), 

U a(l+2) = -ih ■ [R L {d^g{M° cl ,M° c2 })- 1 R R ] a2 , 

V a{i+2) = -m ■ \R R (diag{J^, M°})- 2 R R ] , (A6) 



giving the result 



Ul(i+2) 



Hi a/2 M w cos/3 



M? 



y/2 M W (M* sin/3 + /i cos/3) 

Vi(i+2) = Hi m 



i 



v (|M 2 | 2 + 2il4cos^) 

v 2 (i +2 ) = -turn , (A7) 



where 



M 2 = /x M 2 - M£ sin2/3 = M° M° c2 . 
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These RPV elements correspond to those given in Ref. |7j (with all complex phases ne- 
glected), where the crucial cos/3 dependence of the nonstandard Z°-boson couplings of the 
physical charged leptons = Xi+2) through Ui( i+2 ) is emphasized. The difference between 
U(i + 2)a and U a ( i+2 ) is hence very important. 

The 7x7 Majorana mass matrix for the neutral fermion can be written as 



M 



M 


M z sin^, sin/3 





M 

-M z cosOw sin/? 



-M z smOw cos/3 M z cosO w cos/3 






M z sin6Wsin/? -M z &ax9w cos/3 

—M z cos$^ sin/3 M z cos^ cos/3 

-fio ~Hx -jU a -fis 

-Ho 

-Hi. 

-fa 

-Ha 



(A8) 



with explicit basis (— iB, — iW, W u , ffl d , 1° , 1°, , Z° ). Note that lf d = Z°, and, from the above 
discussion of the charged fermions, we have, for small /ij's, (11,1°, 11) ~ (v e , v»,v^). The 
symmetric, but generally non-Hermitian, matrix can be diagonalized by using unitary matrix 
X such that 



X T M M X = diag{il^o} . 



(A9) 



Again, the first part of the mass eigenvalues, for n = 1-4 here, gives the heavy states, 
i.e., neutralinos. The last part, M^o for n = 5-7 are hence physical neutrino masses at the 
tree level. 
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Table caption 



Table 1 — Numerical 1-loop neutron EDM results from SUSY without R parity, for four 
illustrative cases. All EDM numbers are in ecm. Note that the quark EDM numbers are 
direct output from the numerical program applying our quark dipole formulae; while the 
neutron EDM numbers are from the valence quark model formula, as given in Eq. ([]]). Im- 
parity violating parameters not given are taken as essentially zero. All parameters are taken 
real except those with complex phases explicitly listed in each case, where the real number(s) 
listed then give the magnitude(s). Parameter A here means a common A u and Ad. Only M 2 is 
shown for the gaugino masses; the others are fixed by the unification relationship. Explicitly, 
we use M x = 0.5 M 2 and M s = 3.5 M 2 . The first column under "EDM Results" gives the 
couplings of the loop vertices involved. A g indicates either one of the electroweak gauge 
couplings, while a A' coupling means one with the appropriate admissible flavor indices. In 
the explicit results of the four cases, the latter is always a X 3n . The second column gives the 
reference Feynman diagram figures, when available. The third column indicates whether the 
particular contribution involves a LR squark mixing. In the case that the mixing is involved 
and a R-parity violating (RPV) one is involved in generating a RPV EDM contribution, it 
is marked with "RPV" . 

Figure captions : 

Fig. 1 — Diagram for (i-quark EDM from the gaugino loop. 

Fig. 2 — Diagrams with neutral gaugino-Higgsino mixing for (i-quark EDM. 

Fig. 3 R-parity violating neutralinolike loop diagrams for (i-quark EDM. Naive 
electroweak-state analysis suggests that such a diagram is proportional to a vanishing l° k - 
gaugino mass mixing. 

Fig. 4 — Diagram for (i-quark EDM suggesting involvement of Majorana masses among the 
1° or "neutrinos". 

a 

Fig. 5 — Diagrams for u- and (i-quark EDMs with charged gaugino-Higgsino mixing. 

Fig. 6 — R-parity violating charginolike loop diagram for d-quark EDM. Naive electroweak- 
state analysis suggests that the diagram is proportional to the vanishing l^-W + mass term. 

Fig. 7 — Diagrams for u- and (i-quark EDMs with a fi a mass insertion. 

Fig. 8 — Diagrams for u- and (i-quark EDMs with a B a scalar mass insertion. 

Fig. 9 — Diagrams with a Majorana-like scalar mass insertion for u- and d-quark EDMs. 
For the w-quark case, the direct Majorana-like h u mass insertion is explicitly shown. For 
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the d-quark case, the corresponding direct h d mass insertion is obvious, for a = (3 = 0; 
for a and/or (5 nonzero, the naive direct result from the diagram would vanish, due to the 
vanishing VEVs. 

Fig. 10 — Logarithmic plot of (the magnitude of) the RPV neutron EDM result for 
value between ±2000 GeV, with the other parameters set at the same values as case A in 
Table I. The lines marked by G, C, N, and "Total" give the complete gluino, charginolike, 
neutralinolike, and total (i.e., sum of the three) contributions, respectively. Note that the 
values of the N contributions and those of the C line for < —900 GeV are negative. 

Fig. 11 — Logarithmic plot of (the magnitude of) the neutron EDM result verses tan/3. We 
show here the MSSM result, our general result with the RPV phase only, and the generic 
result with complex phases of both kinds. In particular, the A and /i phases are chosen 
as 7° and 0.1° respectively, for the MSSM line. They are zero for the RPV-only line, with 
which we have a phase of | for X 3U . All the given nonzero values are used for the three 
phases for the generic result (from our complete formulae) marked by "GSSM". Again, the 
other unspecified input parameters are the same as for case A of Table 1. 

Fig. 12 — Logarithmic plot of (the magnitude of) the neutron EDM result versus 9 A (the 
complex phase for the A parameter). The four lines shown are characterized by the tan/3 
values ( 3 or 50 ) used and whether it is for our GSSM result (again with a phase of | for 
X 3n ) — marked by G; or the result for MSSM — marked by M. Again, the A' 3n phase is set 
at | for the G lines, and the /i phase at 0.1° for all; the other unspecified input parameters 
are the same as for case A of Table 1. 

Fig. 13 — Logarithmic plot of (the magnitude of) the neutron EDM result verses M 2 . The 
four lines correspond to the four cases of Fig. 12, each with 9 A set at the dip location, i.e., 
G-3 for GSSM at tan/5 = 3 with 9 A = 2°, M-3 for MSSM at tan/5 = 3 with Q A = -1°, G-50 
for GSSM at tan/3 = 50 with 9 A = 20°, M-50 for MSSM at tan/5 = 50 with 9 A = 3°. All other 
unspecified input parameters are the same as for Fig. 12. 
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Choice of Parameters 
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FIG. 1. Diagram for <i-quark EDM from the gaugino loop. 
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FIG. 2. Diagrams with neutral gaugino-Higgsino mixing for (i-quark EDM. 
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FIG. 3. R-parity violating neutralinolike loop diagrams for <i-quark EDM. Naive elec- 
troweak-state analysis suggests that such a diagram is proportional to a vanishing ^-gaugino mass 
mixing. 
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FIG. 4. Diagram for (i-quark EDM suggesting involvement of Majorana masses among the l° a 
or "neutrinos". 
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FIG. 5. Diagrams for u- and (i-quark EDMs with charged gaugino-Higgsino mixing. 
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FIG. 6. R-parity violating charginolike loop diagram for (i-quark EDM. Naive electroweak-state 
analysis suggests that the diagram is proportional to the vanishing l' k ~W + mass term. 
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FIG. 7. Diagrams for u- and d-quark EDMs with a ^ mass insertion. 
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FIG. 8. Diagrams for u- and <i-quark EDMs with a B a scalar mass insertion. 
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FIG. 9. Diagrams with a Majorana-like scalar mass insertion for u- and <i-quark EDMs. For 
the u-quark case, the direct Majorana-like h u mass insertion is explicitly shown. For the <i-quark 
case, the corresponding direct h d mass insertion is obvious, for a = (3 = 0; for a and/or (3 nonzero, 
the naive direct result from the diagram would vanish, due to the vanishing VEVs. 
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FIG. 10. Logarithmic plot of (the magnitude of) the RPV neutron EDM result for /Jq value 
between ±2000 GeV, with the other parameters set at the same values as case A in Table 1. The 
lines marked by G, C, N, and "Total" give the complete gluino, charginolike, neutralinolike, and 
total (i.e., sum of the three) contributions, respectively. Note that the values of the N contributions 
and those of the C line for ^ < —900 GeV are negative. 
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FIG. 11. Logarithmic plot of (the magnitude of) the neutron EDM result verses tan/?. We 
show here the MSSM result, our general result with the RPV phase only, and the generic result 
with complex phases of both kinds. In particular, the A and phases are chosen as 7° and 0.1° 
respectively, for the MSSM line. They are zero for the RPV-only line, with which we have a phase 
of j for A' 3n . All the given nonzero values are used for the three phases for the generic result (from 
our complete formulae) marked by "GSSM". Again, the other unspecified input parameters are 
the same as for case A of Table 1. 
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FIG. 12. Logarithmic plot of (the magnitude of) the neutron EDM result versus 6 A (the complex 
phase for the A parameter). The four lines shown are characterized by the tan/3 values ( 3 or 50 ) 
used and whether it is for our GSSM result (again with a phase of j for A^) — marked by G; or 
the result for MSSM — marked by M. Again, the X 3n phase is set at j for the G lines, and the /i 
phase at 0.1° for all; the other unspecified input parameters are the same as for case A of Table 1. 
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FIG. 13. Logarithmic plot of (the magnitude of) the neutron EDM result versus M 2 . The four 
lines correspond to the four cases of Fig. 12, each with 9 A set at the dip location, i.e., G-3 for GSSM 
at tan/3 = 3 with 6 A = 2°, M-3 for MSSM at tan/3 = 3 with 6 A = -1°, G-50 for GSSM at tan/3 = 50 
with A = 20°, M-50 for MSSM at tan/3 = 50 with 9 A = 3°. All other unspecified input parameters 
are the same as for Fig. 12. 
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